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Unstable flux flow due to heated electrons in Y-Ba-Cu-O films 
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A flux instability occurs in superconductors at low temperatures, where ee scattering is more 
rapid than ep, whereby the dissipation significantly elevates the electronic temperature while main- 
taining a thermal-like distribution function. The reduction in condensate and rise in resistivity 
produce a non-monotonic current-voltage response. In contrast to the Larkin-Ovchinnikov instabil- 
ity where the vortex shrinks, in this scenario the vortex expands and the quasiparticle population 
rises. Measurements in YiBa2Cu307_^ agree quantitatively with the distinct predictions of this 
mechanism. 
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In a type II superconductor, magnetic fields between 
the lower critical value H c \ and upper critical field H c2 
introduce flux vortices containing a quantum of flux 
$0 = h/2e. Here we have superconducting films in a 
perpendicular applied flux density B, with a transport 
electric current density j in the plane of the film, which 
exerts a Lorentz driving force Fl = j$ - The vortex mo- 
tion generates an electric field E = vB and is opposed by 
a viscous drag i]v (ij is the coefficient of viscosity and v 
the vortex velocity), so that in steady state j$ = rjv and 
the response is Ohmic. Larkin and Ovchinnikov (1986) 
have shown that a dirty superconductor at low tempera- 
tures has a free-flux-flow resistivity related to the normal- 
state value p n by jj] 

p f /p n ~ 0.9B/H c2 (T). (1) 

Approximately the same result, without the precise 0.9 
prefactor, can be obtained by considering the Ohmic dis- 
sipation in the core and temporal changes in the order 
parameter leading to irreversible entropy transfer 0] ; the 
result is also valid for d-wave superconductors that are 
not superclean 0. Eq. [l] is equivalent to r\ m ^qH c2 / ' p n . 
At low levels of j and E, in the assumed dirty limit 
I <C £,Ep/kT c (I is the mean free path and Ep is the 
Fermi energy) , 77 is a constant that is proportional to the 
order parameter A and inversely proportional to the size 
~ £ 2 (where £ is the coherence length) of the vortex |IJ. 
(The Hall effect and transverse component of E are neg- 
ligible for this discussion, as is vortex pinning because of 
the large driving forces 0].) 

At high electric fields and dissipation levels sufficient to 
alter the electronic distribution function and/or the elec- 
tronic temperature, j(E) becomes non-linear and can de- 
velop an unstable region (dj/dE < 0) above some critical 
vortex velocity v*. For the regime near T c , such an in- 
stability has been predicted by Larkin and Ovchinnikov 
(LO) |p[ ; and has been experimentally well established 
[@, g, |§T|lO|- At high temperatures, the electron-phonon 
(ep) scattering time r ep can be shorter than the electron- 
electron (ee) scattering time r ee , preventing internal equi- 



libration of the electronic system and producing a pecu- 
liar non-thermal distribution function. Since the order 
parameter A is especially sensitive to the distribution 
function when close to T c , even moderate values of E 
can sufficiently distort A (via the Eliashberg mechanism 
PI) causing a shrinkage of the vortex core and a removal 
of quasiparticles from its vicinity. This is the gist of the 
LO behavior d, 0. As LO themselves emphasize, the 
effect is most favorable close to T c for superconductors 
with a full gap and, as shown by Bezuglyj and Shklovskij 
@, is dominant for B < B T (with Bt ~ 0.1 T for our 
low-T regime). One of the predictions of the standard 
LO effect is a v* that is B independent. This has been 
confirmed in YiBa2Cu30y_5 in the high-T range || ^0| . 

This work investigates the opposite regime of T <C T c 
and B > Bt, where A is not sensitive to small changes 
in the distribution function. Furthermore because r ee < 
T ep as T — > 0, the distribution function remains thermal 
like and the electronic system suffers mainly a tempera- 
ture shift with respect to the lattice §, |S| @. Then 
instead of the standard LO picture described earlier, we 
consider a more transparent scenario where the main ef- 
fect of the dissipation is to raise the electronic temper- 
ature, create additional quasiparticles, and diminish A. 
The vortex expands rather than shrinks, and the viscous 
drag is reduced because of a softening of gradients of 
the vortex profile rather than a removal of quasiparti- 
cles. This sequence of events is almost opposite to the 
standard LO picture and represents a new type of unsta- 
ble regime prevalent at T <C T c . All experimental mea- 
surables can be calculated without ambiguity, and the 
predicted field dependencies and full j(E) curves fit the 
experimental results without any adjustable parameters. 

Previously some deviations from LO behavior at inter- 
mediate temperatures — such as a -B-dependent v* — were 
treated through modifications to the LO effect, such as an 
intervortex spacing 1^ that exceeds the energy-relaxation 
length l e jl^] or by inclusion of thermal effects |l2|, 
Those treatments do not apply to the present T <C T c 
regime where l e ~ yjDr ep ~ 100-1000 nm is larger than 
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Z = 1.075-v/CTB ~ 10 " 50 nm (-D = 3 x ICT 4 m 2 /s is 
the diffusion constant 

Here we take a "bottom-up" approach and start from 
the T ~ limit: The total input power jE travels from 
electrons to lattice and from there to the bath, so that 
To < T p < T' , where To an d T p are the bath and phonon 
temperatures, and T' is the raised non-equilibrium elec- 
tronic temperature. Macroscopic heating, represented by 
T p - T = RthjE, is < 5% of the total increase T- T for 
the worst case dissipation so that T p « T (Here i?th~l 
nK.cm 3 /W is the total thermal resistance between the 
film and the bath; see experimental section. It will be 
seen later that the specific-heat integral heavily weights 
the higher temperatures, so that T p - Tq is quite negligi- 
ble.). The energy relaxation between electrons and lat- 
tice occurs by inelastic ep scattering, and is characterized 
by an effective time r e ^<r ep >. The principle contribu- 
tions to p n come from impurities and phonons. Since the 
phonon temperature remains near the bath, /o n 's value 
will not change as the non-equilibrium T' rises. Thus 
putting p n (To) and H C 2(T') (since H c2 does depend on 
T') into previous equations gives the j(E) response in 
terms of T': 

j = vr)(T')/$ = EH c2 (T')/p n (T )B. (2) 

One now has the ingredients for calculating the critical 
field dependencies in a few steps. The j(E) function of 
Eq. |^ is non-monotonic since E is multiplied by H C 2(T') 
(or rj(T')) which drops rapidly to zero as T'— >T C with in- 
creasing dissipation. The instability occurs at dj / dE = 0, 
which happens at a certain value T'—T* where H C 2(T') 
(or rj{T')) drops sufficiently rapidly. T 1 depends explic- 
itly only on the power density jE = nv 2 B/$o and on 
quantities that depend on T' itself (r e , specific heat, etc.). 
Hence at the instability, j* E*=v* 2 r/(T*)B /$o= constant, 
which gives the critical-parameter field dependencies: 

v* cx l/VB,E* cx y/B,f cx 1/Vb, and p* cx B (3) 

using E*=v*B, j*=v*rj(T*)/$ , and p*=E*/j*. This 
gives v* cx 1/ \f~B in a natural way, consistent with our 
measurements in this regime and in contrast to the B- 
independent v* of the pure LO effect near T c . 

To derive the complete j(E) response and absolute val- 
ues of critical parameters we calculate T' and insert it 
into Eq. 0. We take iT c2 (0)=120 T (T|] with the WHH 
function |17tl for interpolation between this iJ C 2(0)=I20 
and -ff C 2(T c )=0 (There is some theoretical controversy re- 
garding the exact form of H C 2(T); however, empirically, 
direct measurements of H C 2{T), within their uncer- 
tainly, seem not to depart drastically from the WHH 
function, and the exact functional shape does not cru- 
cially affect our conclusions.). The dissipation raises the 



electronic energy by jEr e , which is related to T' by 

T' T' 

jEr e « E2Hc ^ T 2 Te « / < T ) dT « / < T ) dT ( 4 ) 

Pn{T )B J J 

where c(T) is the electronic specific heat. To calcu- 
late c(T), Y!Ba2Cu 3 07_5 is modeled as a layered d- 
wave superconductor: A^(T) = A (T)[fc 2 — ky]/k 2 ps 
A (T) cos(26'), taking the BCS temperature dependence 
for Ao(T) and Ao(0)=19 meV from tunneling and in- 
frared measurements Q. Then c(T)=d{J2k E kfk}/dT, 
where fk—[exp(Ek/ksT) + 1] _1 is the Fermi-Dirac dis- 
tribution, E k = VCfe +A?., Cfc = £fe - M) and ^=0-2 eV is 
the chemical potential [|I9|. With the replacement of ^2 k 
by J d(N(0) J d9/2ir (where ^(0) is the normal density 
of states at /z) this c(T) is now inserted into Eq. ||, which 
is integrated numerically to obtain T' and thence j(E) 
from Eq. |[ 

The numerical results of the above calculation are 
shown in Fig. |l|. The j(E) curves of panel (a) have 
an infinite slope at the instability (j*, E*), indicated by 
arrows, and then exhibit negative slope. This negative 
sloped portion is experimentally forbidden in a current 
biased measurement and instead will be manifested as a 
vertical jump in E. The electronic temperature T' rises 
from the bath value To to T* at the instability. Panel 
(b) shows the computed T*, which is independent of B 
as expected but has a slight dependence on the bath tem- 
perature Tq. 
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FIG. 1: Numerical results obtained from solving Eq. W, as 
described in the text, (a) Theoretical j(E) curves, at To= 
OK. The onset of negative slope, indicated by arrows, marks 
the instabilities, (b) The critical temperature T* for different 
initial temperatures To; T* is independent of B and r e . 

In order to conveniently scale the experimental curves, 
the exact numerical j(E) function derived above — and 
plotted in Fig. |l| — can be cast into a mathematically 
more manageable form by noting that the R.H.S. of Eq. |] 
can be approximated [? ] by jEr e w A(T')n q (T') — 
A(T p )n q (T p ) « A(T')n q (T'). A and n q are connected 
through the statistical equations of the previous para- 
graph. Numerically computing the number of quasiparti- 
cles excited above the gap (taking the anisotropic d-wave 
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gap together with a BCS temperature dependence as dis- 
cussed above) we obtain the following d-wave generaliza- 
tion of the A — n q relationship [? ]: (A/A ) 2 = f(n q /n), 
where n — 2.7 x 10 21 cm~ 3 is the carrier concentration 
and f(x) ~ 1 - 0.4386a; - 1539x 2 -I- 40381a; 3 - 345217a; 4 
(despite the appearance of successively increasing co- 
efficients in /(a;), the terms rapidly converge because 
x=riq/n~ kBT/Ep ~ 0.01). Combining this with the 
earlier jEr e « A(T')n q (T'), v = j<t> /v = j$ B/E, 
■q rs H C 2$>o/ p n , and 77 cx A, we get a convenient closed 
form for the non-linear j(E) characteristic: 

with x = n q /n = 0.0245 x E 2 /E* 2 , and E* = 
v /0.0245p„(To)BnAo/r £J ff C 2(T ) is the value of E at the 
instability peak. We now turn to the experimental details 
and results. 

The samples are c-axis oriented epitaxial films of 
YiBa 2 Cu 3 07_5 on (100) LaA10 3 substrates with T c 's 
around 90K and of thickness t f» 90 nm. Electron-beam 
lithography was used to pattern bridges of widths 10 f=a 
2-20 /im and lengths I » 30-200 (tm. Altogether ten 
samples were studied at 12 temperatures (1.6, 2.2, 6, 7, 
8, 10, 20, 27, 35, 42, 50, 80 K) and at 13 flux densities 
(0.1, 0.2, 0.5, 1, 1.5, 2, 10, 11, 13, 13.5, 13.8, 14, 15.8 T). 
The electrical transport measurements were made with 
a pulsed constant current source, preamplifier circuitry, 
and a digital storage oscilloscope. The pulse rise times 
are about 100 ns with a duty cycle of about 1 ppm, which 
for the narrowest bridges result in effective thermal re- 
sistances of order 1 nK.cm 3 /W. Note that the j values 
in the experiment are an order of magnitude lower than 
the depairing current density p0[ and the applied flux 
densities exceed the self field of the current by at least 
two orders of magnitude. Further details about the ex- 
perimental techniques are discussed elsewhere pl| . 

Fig. ||(a) shows a typical set of experimental j(E) 
curves. The last stable datapoint (j* , E*) of each curve 
is at the tail of each arrow. The slightest further increase 
of j >j* causes a drastic vertical jump in E as shown 
by the arrows (the voltage pulse jumps off the scale of 
the oscilloscope upto the compliance limit of the current 
source). The jumps show only a small hysteresis < 3% of 
j*. (As expected for a current-biased measurement, the 
break occurs slightly before the slope has become quite 
vertical |fi2[.) Fig. 0(b) shows the same data plotted as 
the R.H.S. of Eq. Qvs j — j c o, the excess current den- 
sity over the depinning value. The data scale well and 
tend toward homogeneous linearity. Note that the col- 
lapse implies an excellent proportionality between p and 
B over the entire range. 

Fig. ^ shows experimentally measured B dependencies 
ofw* and p* for To=1.6 K, demonstrating excellent agree- 
ment with Eq. H (the other dependencies E* cx \[~B and 




FIG. 2: Experimental current-voltage characteristics, (a) 
Raw values of E versus j at T=20K and applied B values 
of (from lowest to highest curve) 3, 5, 8, 11, 13.8, and 15.8 
T. The last symbol on each curve is right at the instabil- 
ity. Slightest further increase of j (entering the forbidden 
negative-sloped portion of the theoretical curves of Fig. 1) 
causes E to make discontinous vertical jumps (arrows), (b) 
The same data plotted as E + = (E/B) x y/ /(0.0245.E 2 /E* 2 ) 
(as per Eq. ^) versus j —jco, where the critical depinning cur- 
rent density j c o is defined at E + — 1 V/cm.T (the scaling is 
not affected by the choice of criterion). 




FIG. 3: Variation of critical parameters with flux density. 
The measurements were made at T=1.6 K. (a) The critical 
velocity shows a v* oc 1/yB trend (solid line is a 1/vB fit), 
(b) The critical resistivity is proportional to the flux density 
(straight line is a guide to the eye). 

j* cx 1/VB follow from p* cx B and v* cx 1/VB). The 
v* cx \j\J~B dependence was found to be ubiquitous for 
all of our low-T measurements in ten samples (spanning 
1.6 K < T < 50 K and 0.5 T < B < 15.8 T) and has 
also been seen by Xiao et al. |l(J at intermediate tem- 
peratures (at the lower end of their T~60-90 K range 
of study). Note that the excellent linearity between be- 
tween p* and B demonstrates the independence of i] on B 
in this regime; then the resistivity is simply proportional 
to the number of vortices and hence B. 

The final step in the analysis is to extract r e from 
the data. Inverting Eq. ^ evaluated at the instability 

(i.e., at T'=T*), we get r e « h p ^]e^ / c(T)dT; the 

To 

T*(T ) function comes from the model itself (Fig. 0(b)). 
T e calculated in this way is shown by the circles in 
Fig. ||. If the data are analyzed in the LO frame- 
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work (r e « D[14C(3)] 1 / 2 B 2 /( 7r - E '* 2 )) a markedly longer 
r e is obtained shown by the squares. For compar- 
ison the arrow shows the theoretical T— ► limiting 
value estimated from spontaneous emission of phonons. 
(Emission of non-equilibrium phonons dominates ep scat- 
tering by thermal phonons below a cross-over tem- 
perature T x = mvFC s /kB~43 K producing an essen- 
tially T-independent scattering rate (since kT<^Ep) of 
T ep ={3Trp M Csh 2 )/(2C 2 mk F ) » 8.4 x 1CT 11 , where C*«6x 
1CK 13 ergs is the coupling between electron energy and 
crystal dilation, c s — 4 km/s is the sound velocity, and 
Pm=6 g/cm 3 is the mass density [22 ; for further discus- 
sion see, for example, Kittel pq].) 
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FIG. 4: The energy relaxation time r e extracted from the 
measured data using the present low-temperature instability 
model (circles) and from LO theory (squares). The leftward 
arrow represents the theoretical T — > limiting value of r ep 
from phonon emission (please see text). 

In conclusion, we investigated the low-temperature 
regime of flux motion driven far beyond free flux flow, 
and observed an instability under all conditions of fields, 
and temperatures from ^T c /2 down to essentially T ~ 0. 
The nature of this low-temperature instability seems to 
be well described by a model where the electron gas is 
heated above the phonon temperature leading to the gen- 
eration of quasiparticles and loss in viscosity as the vortex 
core expands and A is reduced. This scenario is differ- 
ent from the standard LO picture (dominant mainly near 
T c ) where the vortex shrinks and quasiparticles leave its 
vicinity. Because the present effect prevails even at tem- 
peratures well below T c (where most superconductive de- 
vices operate) it becomes an important consideration in 
the design of applications where the superconductor op- 
erates in the dissipative regime ^4|, since the instability 
triggers an abrupt rise in dissipation at j values much 
below the depairing current density. Detailed predic- 
tions of the model, including field dependencies of critical 
parameters and shapes of j(E) curves were experimen- 
tally verified, and the effect provides an estimate of the 
T^O time scale for energy exchange between quasiparti- 
cles and phonons. 
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